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Abstract: Observations of the Crab Nebula have proven to be the best tool to calibrate and to characterize 
the performance of a Cherenkov telescope. Scientifically, it is interesting to measure the energy spectrum 
of the Crab Nebula close to the inverse-Compton peak where a deviation is expected from the power 
law seen at energies above 300 GeV. Additionally, it is important to search for pulsed emission from the 
Crab Pulsar at energies beyond the 10 GeV upper limit of the EGRET pulsar detection. Since current 
models predict a cut-off in pulsed emission between 10 and 100 GeV, measurements at energies close 
to this range may help to discriminate between them. We observed the Crab extensively in the 2006- 
2007 season during the VERITAS 2- and 3-telescope commissioning phases. Using this data set we 
reconstructed a preliminary energy spectrum of the signal from the Crab Nebula. We also measured the 
optical pulsed signal to validate our GPS time-stamping and barycentering techniques and obtained an 
upper limit for the pulsed emission at gamma-ray energies. 



Introduction 

VERITAS, located in southern Arizona (32°N, 
111°W) at an elevation of 1268 m, is an obser- 
vatory designed to explore the very high energy 
(VHE) gamma-ray sky in the energy band between 
100 GeV and 50 TeV [13]. It consists of four 12m 
diameter Cherenkov telescopes, each with a 499- 
pixel camera. It uses the stereoscopic imaging at- 
mospheric Cherenkov technique (lACT) to detect 
gamma rays from astronomical sources. VERITAS 
has moved through phases of 2-, 3- and 4-telescope 
operation during its construction and commission- 
ing this past year. We observed the Crab Neb- 
ula and Pulsar extensively during the 2- and 3- 
telescope phases. This paper discusses the anal- 
ysis of these data. We also present a preliminary 
differential energy spectrum of the nebula and a 
preliminary upper limit for the pulsed gamma-ray 
emission. 

The Crab Nebula and Pulsar 

The Crab Nebula is a plerion-type supernova rem- 
nant (SNR). It was the first source definitively 



detected at TeV energies, by the Whipple col- 
laboration [16]. Its spectral energy distribution 
is well measured across a broad range of wave- 
lengths and it is considered the standard candle in 
TeV astronomy due to its strong, stable emission. 
The nebula emission is believed to be powered by 
the conversion of the rotational slowdown of the 
central pulsar to synchrotron radiation. Inverse- 
Compton scattering of those synchrotron photons 
with the primary accelerated electrons gives rise 
to the emission in the gamma-ray regime as ex- 
plained by the Synchrotron-Self-Compton model 
[6]. The Crab Nebula spectrum has been measured 
independently by many groups above ^--^200 GeV, 
with good agreement [9, 1, 15, 2]. Consistency 
with this established spectrum is an important test 
for any new gamma-ray telescope. 

At the center of the Crab Nebula lies the Crab Pul- 
sar (PSR 0531 H-21). Periodic emission from the 
pulsar, with a 33 ms period, has been observed 
in all energy regimes with the exception of very 
high energy gamma rays. The EGRET detection of 
pulsed emission in the 100 MeV - 10 GeV energy 
band, with no evident turnover of its spectrum, is 
currently the highest energy detection [14]. None 
of the ground-based gamma-ray experiments sen- 
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sitive in the adjacent VHE domain has yet provided 
any detections. Currently there are two primary 
models attributing the pulsed emission to curvature 
and synchrotron radiation from the energetic elec- 
trons/positrons accelerated in the strong magnetic 
fields of the pulsar magnetosphere. These mod- 
els differ by whether the emission region is near 
to the magnetic poles of the neutron star (the polar 
cap model [5]), or near to the light-cylinder at the 
outer magnetosphere (the outer gap model [10]). 
In the polar cap model, a super-exponential cut- 
off occurs in the pulsed spectrum at a few GeV 
due to the magnetic pair production in the strong 
field close to the pole. On the other hand, in the 
outer-gap model a cut-off at a few tens of GeV 
is expected due to photon-photon pair production, 
which has a weaker energy dependence. By detect- 
ing the pulsed emission or placing stringent upper 
hmits in the GeV-TeV energy range, one can po- 
tentially discriminate between these models. VER- 
ITAS has a better sensitivity at lower energies than 
previous-generation lACT instruments. This sug- 
gests the possibility of detecting the pulsed emis- 
sion, if an outer gap model applies, or of strongly 
constraining both models if pulsed emission is not 
detected. 

Data Acquisition and Analysis 

The data set used in this analysis was selected from 
the observations made using 2 and 3 telescopes be- 
tween October 2006 and February 2007. During 
that time, different trigger and observing configu- 
rations were tested to determine the optimum con- 
figuration for the array. Only data with this op- 
timum configuration are used in the present anal- 
ysis. Data-quahty selection, based on the stabil- 
ity of trigger rates, response of the instrument to 
background cosmic rays and the deadtime of the 
system, was additionally used to reject low-quality 
data. The final data set, after this quality selection, 
is 15.7 hours of 2-telescope data and 3.3 hours of 3- 
telescope data. Data on both sets were collected in 
the so-called "wobble" observing mode, for which 
the telescopes are pointed such that the source is 
at a 0.5°offset from the position of the camera cen- 
ter. For 2-telescope data we required that both tele- 
scopes trigger within a certain coincidence time. 



while for 3-telescope data we required that any two 
out of three telescopes trigger to record the event. 

The data were analyzed using the standard 
VERITAS analysis package [4, 3] in which cali- 
bration, cleaning and parametrization of the signal, 
reconstruction of the shower parameters and the 
gamma/hadron separation cuts are applied to the 
data. The quality-selection cuts which require the 
total recorded signal (Size) to be greater than 400 
d.c, the number of pixels in the image (^o fPix) 
to be greater than 5 and the distance of its centroid 
from the center of the camera (Dist.) to be in the 
range (0.05°, 1.3°) are applied before the stereo 
reconstruction. The gamma/hadron separation cuts 
based on the mean scaled width and length (0.05 < 
MSW < 0.95 and 0.05 < MSL < 1.15 re- 
spectively) and the squared distance between the 
source and the reconstructed gamma ray position 
(0'^ < 0.025) are applied at the final step. The 
same data were also analyzed using an independent 
analysis package to confirm the results. 

Detection 

The Crab Nebula is strongly detected during both 
2-telescope and 3-telescope phases. The detec- 
tion is at the 7 7/min level with a significance of 
31 a/^/hr for three telescopes and at 3.4 7/min 
level with a significance of 20 a/ Vhr for two tele- 
scopes using the cuts described above. The ex- 
cess 9'^ distributions for each data set can be de- 
scribed by a sum of two Gaussian functions [2]. 
The fit to 3-telescope data gives Aaei = 0.17 with 
cri = 0.067 and (72 = 0.13. 

Energy Spectrum 

The energy for each succesfully reconstructed 
event is estimated using lookup tables based on 
the Size, the impact parameter of the event and 
the zenith angle. The lookup tables are generated 
using the simulations of the gamma-ray showers 
and the detailed detector response. Effective col- 
lection area of the instrument as a function of en- 
ergy of the event, the zenith angle and offset of 
the source from the pointing direction is also de- 
termined using the same simulations. To calcu- 
late the differential flux from the Crab Nebula, a 
looser gannma/hadron seperation cuts are apphed 
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to the data to keep more events in lower energy 
bins. In each energy bin, the gamma-ray events 
passing the cuts are weighted by the inverse of the 
effective area of the instrument for the energy of 
the event, and the difference of this weighted sum 
of the gamma-ray events from the on source region 
and the normalized weighted sum of the gamma- 
ray events from the background regions is calcu- 
lated. The differential flux for each energy bin is 
determined by normalizing this difference by the 
live-time for that energy bin and the bin width. 
Figure 1 shows the measured energy spectrum of 
the Crab Nebula using the 3-teIescope data set fit- 
ted to a simple power law between 250 GeV and 
6 TeV. The minimum energy bin in the fit is de- 
fined by requiring the energy bias lower than 10% 
and the maximum energy bin is defined by requir- 
ing a detection of minimum of 2 sigma in each bin. 
The spectrum is checked for consistency using two 
independent analysis packages and it also agrees 
very well with results from other experiments. 
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Prob 0.2102 

Norm 3.626e-07± 1.516-08 

Index -2.544 ±0.04918 
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Figure 1: The reconstructed Crab Energy Spec- 
trum using 3-telescope data set. 



of the first radio pulse at the solar system barycen- 
ter The radio ephemerides for the Crab pulsar are 
obtained from the Jodrell Bank Observatory [12]. 
Barycentering libraries used in this analysis were 
developed [7] and tested for this analysis using the 
optical data taken with a special set-up explained 
in the following section. The same data set is also 
tested for periodicity by an independent temporal 
analysis [11]. 

The Crab Optical Pulsar 

Before analyzing the gamma-ray data for period- 
icity, we took special runs using only the central 
pixel of one telescope to detect the optical pulsed 
signal from the Crab pulsar. This detection of 
the periodic optical signal confirms that our on- 
line time stamping and our offline data analysis are 
capable of detecting the 33 ms pulsed emission. 
For the special optical runs, we ran the system at 
the single photo-electron trigger level, triggering 
the system at 20 kHz and recording only the time- 
stamps and number of pulses collected between 
each trigger from the central pixel and a pixel at 
the outer edge of the camera. Figure 2 shows the 
phase plot of the signal obtained by combining the 
data from 3 runs taken on 2 different nights, for a 
total of 30 minutes. It clearly shows the main and 
inter pulses of the optical pulsed emission. 
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Temporal Analysis 

The temporal analysis of the Crab signal is used 
to search for pulsed emission from the Crab Pul- 
sar. The event arrival times are tagged online by 
GPS clocks at each telescope and are transformed 
to the solar system barycenter to compensate for 
the motion of Earth. We then calculate the relative 
phase of the pulsar rotation for each event with re- 
spect to a reference epoch which is the arrival time 



Figure 2: The phase histogram of the optical signal 
detected by VERITAS. 



The Pulsed Emission Upper Limit 

Since the pulsed emission is expected only at the 
lowest energies to which VERITAS is sensitive, the 
Size-cut used in the standard VERITAS analysis 
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removes all relevant events. Thus, to search for 
the pulsed emission, all data which pass the stan- 
dard quality cuts with the exception of the Size- 
cut are used. Since the angular resolution is de- 
graded at low energy, the standard 6*^ -cut was also 
relaxed to 9^ < 0.04 to enhance the acceptance of 
low-energy events. The binning of events by phase 
acts as a separation of signal and background, and 
no further gamma/hadron separation cuts are ap- 
plied. These loose cuts increased the overall num- 
ber of candidate gamma-ray events by a factor of 
27. The significance of the pulsed excess is calcu- 
lated by assuming the pulse shape as seen in the 
EGRET energy range, with a defined pulsed emis- 
sion region between phases [—0.06 — 0.04] and 
[0.32 — 0.43] for the main pulse and inter-pulse re- 
gions respectively. The remaining phase regions 
are used to estimate the background for the pulsed 
emission. No evidence for pulsed emission is seen. 
An upper limit for the pulsed emission is calculated 
using the method of Helene [8], assuming Gaus- 
sian statistics. The upper limit on the integral flux 
of the pulsed signal is then derived separately for 
the 3-telescope and 2-telescope data. The integral 
flux upper limit obtained using 3-telescope data set 
is 8.17 X lO^^^cm^^s^^ above the energy thresh- 
old of 200 GeV and the upper limit for 2-telescope 
data set is 3.5 x 10~^^cm~^s~^ above 240 GeV as 
plotted on figure 3. 
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Figure 3: The integral flux upper limits for the 
pulsed emission obtained from 3-telescope and 2- 
telescope data set are shown with filled circles 
along with the results from other experiments. 



Conclusion 

Steady emission from the Crab Nebula is detected 
at very high significance by VERITAS. The steady 
signal from the Crab was used to perform vari- 
ous performance checks on the system for differ- 
ent configurations. A preliminary differential spec- 
trum of the Crab was measured. A simple power 
law fit to the data points gives ^ = (3.63 ± 
0.15) X 10-7(^)(-2-54±o.O5) TeV-^m-\s-\ 
Pulsed emission from the Crab Pulsar is also 
searched for with the same data set, using looser 
cuts. Integral flux upper limit for the pulsed emis- 
sion of 8.17 X lO^^^cTO^^s^^ above the energy 
threshold of 200 GeV is calculated from the 3- 
telescope data set. 
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